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Abstract 

We investigated //"'" decays at rest produced at the ISIS beam stop target. Lepton flavor (LF) 
conservation has been tested by searching for Vf, via the detection reaction p ( ) e"*" ) n . No i>e. 
signal from LF violating /x+ decays was identified. We extract upper limits of the branching 
ratio for the LF violating decay /x"*" — > e"*" + i>e + compared to the Standard Model (SM) 
/i+ ^ e+ + i^e + decay: BR < 0.9(1.7) • lO^^ (90%C.L.) depending on the spectral distribution 
of characterized by the Michel parameter p = 0.75(0.0) . These results improve earlier limits 
by one order of magnitude and restrict extensions of the SM in which i>e emission from decay 
is allowed with considerable strength. The decay p,~^ — > e^ + z^e + i^^t as source for the z^e signal 
observed in the LSND experiment can be excluded. 

PACS numbers: 13.35.Bv, 14.60.St, ll.30.Fs 



INTRODUCTION 



In the Standard Model (SM), the main decay mode of positive muons is the decay into 
a positron and two neutrinos fi^ e+ + z/ + z/'. Assuming conservation of the additive 
lepton family or flavor (LF) numbers L^, and L^, the neutrino flavors are fixed to be z/ = 
and u' = v^. The neutrinos are massless with the Vf, being a left-handed neutrino, the a 
right-handed anti-neutrino. The structure of the muon decay can be described by the V-A 
theory of weak interactions. Therefore /x decay as a purely leptonic process has been used to 
study with high precision the SM of weak interactions. The Lorentzian V-A structure of the 
u"*" decay can be tested by measuring the massive leptons, i.e. the initial ii^ and the final 

n □ 

e"*" m or by investigating the neutrino energy spectrum However, to test conservation 
of the LF numbers and in ^ decay it is essential to observe the final neutrino states 
P]. All tests so far show no deviations from the SM. 

However, the LF number violating decay mode /i"*" — > e"*" -|- ^ -f is allowed in many 
extensions of the SM, e.g. left-right (LR) symmetric models P, QI^Qi GUT models with 
dileptonic gauge bosons extensions involving additional scalar multiplets or super- 
symmetric models with R parity violation jl^, together with the LF number violating decay 
/i+ — i> e"*" + 7 3| • Although the energy scale of LR symmetry of weak interactions or the 
appearance of supersymmetric particles is expected to be in the range of 0.1-1 TeV, precision 
measurements at intermediate energies can provide essential restraints on the parameters 
used in various models. Therefore, the detailed investigation of the /i"*" decay plays a major 
role in determining the structure of weak interactions and the precision of lepton number 
conservation. 

On the other hand, there are clear evidences for neutrino oscillations from experiments on 
atmospheric, solar and reactor neutrinos Since v oscillations violate the conservation 
of the lepton family numbers, such results enhance the interest of searching for direct LF 
number violation. In addition, there is a positive Ve signal from the accelerator experiment 
LSND Ql which could be explained a priori as an indication for z/^— oscillations of z/^ from 
jJi^ ^ e"*" + z/g + z/^ or directly for the decay mode /i"*" — > e"*" + z/g + ^z/^. Due to limited 
statistics and energy resolution, this ambiguity is not resolved by the LSND experiment 
itself. 

The spallation source ISIS at the Rutherford Laboratory, UK, is a unique source of /x"*" 
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to study such decays. The KArlsruhe Rutherford Medium Energy Neutrino experiment 
investigated the neutrinos produced at ISIS through the decays of 7r+ and yU+ at rest. One 
purpose of the KARMEN experiment was the investigation of //-nucleus interactions on 
^^C [15]. The good agreement of the measurements with theoretical predictions allowed a 
sensitive search for processes forbidden in the SM such as z/-oscillations, z/^— Ve and v^,-^ v^. 
in the appearance mode |3] and z/g— i/^. in the disappearance mode |13| or non-SM decay 
modes of vr"*" and 

In this letter we report the results of the search for v^. from /i"*" decay at rest (DAR). 
Note that z/e from non SM interactions can be produced by either ^ e"*" + z/e + ^z/-* 
in the ISIS target or by oscillations Vp,-^ Ve of on their way to the detector with z/^ 
being produced at ISIS in SM /i"*" decays. While the z/g energy spectrum is fixed in the 
DAR /i^ — * e"*" + z/g + with a spatial flux according to a dependence, the energy 
and spatial distributions of z/g from oscillations strongly depend on the oscillation param- 
eters, i.e. the mass difference Am^^- = — m||. With its excellent energy resolution 
of aE/E = ll%/y^£'[MeV], the KARMEN detector is able to separate different scenarios 
for potential z/g occurrence. Although the search for oscillations z/^— ^ z/g and for the decay 
/i+ ^ e^ + z/g + '^z/^ use the same detection reaction p ( z/g , e"^ ) n for z/g, the different physics 
and consequently the different e"*" spectral distributions result in two separate analyses. 



z^g FROM LF VIOLATING /i+ DECAYS 

In the SM, applying the V-A theory, the energy spectra of massless neutrinos from /i"*" 
decay //^ — e"^ + z/g + z/^ can be calculated neglecting radiative corrections as [3] 

N(e)de oc [ 3(1 - e) + -p(4e - 3) ]de (1) 

3 

with the relative energy e = E^^/Emax, Emax = 52.83 MeV for the decay at rest, and the 
Michel parameter p = 0(0.75) for Ve{i^p), respectively. 

Looking for physics beyond the SM, z/g's from /i^ — > e+ + z/g + *^z/^ in general have non- 
zero mass and contributions of left- or right-handed chirality eigenstates. As only the 
z/g in the decay p^ — > e^ -|- z/g -|- '-z/-' is identified in the experiment, the second emitted 
(anti)neutrino '•z/^ is not determined. Since our experimental result sets upper limits on 
e"*" + z/g -|- these limits also apply for the specific case '•z/^ = z/^, which is the 
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dominant one for certain model assumptions [? 1. 

n 

Taking the actual direct mass limits for z/g, ra{ve) < 2.2 eV [18[ (and hence for z/^ and Vr 
masses through the mixing manifested in the experiments on neutrino oscillations), the v 
masses are very small compared to the mass of the charged leptons or the energy scale of the 
neutrinos from /i+ decays at rest. Assuming Majorana type neutrinos, the v^, of left-handed 
chirality emitted in the SM decay /U"*" — > e"*" + i/g + t'/x could be detected via the reaction 
p ( z/g , e"*" ) n since there is no distinction between Ve and Ve ■ However, the detection of 
emitted in muon decays with left-handed helicity would be strongly helicity-suppressed 
(1— /5 ~ o(10^^^) for a neutrino of 10 MeV energy and a rest mass of 2 eV/c^) since only right- 
handed anti-neutrinos are absorbed via p ( z/g , e"*" ) n . In the case of Dirac type neutrinos, 
the above argument applies for left-handed chirality states of the z/g emitted. Therefore, 
KARMEN as any detector up to date is sensitive only to right-handed z/g's. 

With the rest masses much smaller than the energy of all neutrinos emitted in 
/i"^ — i> e"*" + z/e + *^z/\ an analytical description of the neutrino spectra similar to the one 
in equ. can be applied, with the spectral parameter p to be specified, replacing the SM 
Michel parameter p. In some SM extensions with p^ ^ e^ + ^ + z/^, the Ve and take 
the places of the SM and z^e, respectively, with p(z/e) = 0.75 jsf. In others, p = for the 
emitted z^e j^- In our analysis, we therefore investigate the Ve emission for a variety of p 
parameters. 



EXPERIMENTAL CONFIGURATION AND DATA EVALUATION 

The experiment was performed at the neutrino source of the ISIS synchrotron which 
accelerates protons to an energy of 800 MeV before striking a massive beam stop target. 
On average, 4.59 ■ 10~^ vr"'' per incident proton are produced which are stopped within the 
target and decay at rest. Neutrinos emerge isotropically from the consecutive decays at rest 
(DAR) 7r+^ /^^+ and p^ ^ e"*" + z/e + 20| assuming the z/-fiavors of the SM decay 
channels. Neutrinos from p~^ DAR have a continuous energy spectrum according to equ. (P). 
Due to the narrow time structure of 525 ns of the proton pulses muons are produced in a 
short time window compared to their lifetime of 2.2 ps. 

The neutrinos are detected in a 56 1 scintillation calorimeter [21] at a mean distance of 
17.6 m from the ISIS target. The calorimeter is a mineral oil based scintillator segmented 
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into 512 independent modules. Gadolinium within the module walls allows effective neutron 
detection via Gd ( n,7 ) in addition to the capture on the hydrogen of the scintillator via 
p(n,7). The scintillation detector provides an almost pure target of ^^C and for v- 
interactions. Three veto layers ensure a search for LF violating yU+ decays almost free of 
cosmic background. 

i>e's from /i+ decay can be detected via the (e^, n) sequence from charged current reactions 
p(z/e,e+)n and ^^C ( z/g , e"*" n ) ^^B. Hence, the signature is a prompt e"*" and a delayed, 
spatially correlated 7 signal from the capture of the thermalized neutron by p ( n,7 ) with 
= 2.2 MeV or Gd(n,7) with E ^7 = 8MeV. The flux averaged (taking equ. (0) with 
p = 0.75) cross section of p ( z/g , e"*" ) n is a = 93.5 ■ 10^'*^ cm^ j22]. The ""^^C ( Uej^'n. ) ^""^B 
contribution to {e~^,n) sequences has a cross section of a = 8.52 ■ 10"^^ cm^ j23| which 
is further reduced relative to p(z/e,e^)n by the abundance ratio H/C=1.767 within the 
scintillator. 

A positron candidate is accepted only if there is no activity in the central detector or in 
the veto system up to 24 /is beforehand. The prompt event is searched for in an interval of 
0.6 fis to 10.6 fis after beam-on-target. The time structure of the prompt e"*" event relative 
to the proton pulses has to follow the fi^ decay time constant of 2.2 fis. The expected 
visible e"*" energy has been simulated in detail based on Ue spectra with different values of 
the parameter p including both detection reactions p ( i/g , e"*" ) n and ^^C ( z/g , e"^ n ) ^^B. As 
a result, the prompt energy is required to be within 16 MeV < E(prompt) < 50 MeV (see 
Fig.P). 

The time difference between the e^ and the capture 7 is given by the thermalization and 
capture of neutrons and can be approximated by an exponential with a time constant of 
Tn ~ 120 /is. Therefore, the delayed event has to appear within 1.3 m^ around the prompt 
event position, correlated in time (5 /is < At < 300 /is) with a visible energy E(delayed)< 
8 MeV. For further details of the data reduction of sequential event signatures and of the 
neutron detection in KARMEN see also Q]. 

The raw data investigated in this search were recorded in the measuring period of Febru- 
ary 1997 to March 2001 and represent the entire KARMEN 2 data set which corresponds to 
9425 C protons on target with 2.7 ■ 10^^ decays in the ISIS target. Applying all evalua- 
tion cuts, 15 candidate sequences remain with prompt energies as shown in Figure Q The 
expected background amounts to 15.8 ± 0.5 events. This number comprises 3.9±0.2 events 
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FIG. 1: Visible energy distribution of candidate events with background expectation (shaded area). 
The sohd and dashed hues show the 90%C.L. Umit for an additional i>e. signal with a spectral 
parameter p = and p = 0.75, respectively. 

from cosmic induced sequences as well as v induced reactions such as intrinsic source con- 
tamination of z/g (2.0±0.2), Ve induced random coincidences (4.8±0.3) and (e~,e+) sequences 
from ^^C ( z/g , e~ ) ^^Ng.s, with subsequent ""^^N decay (5.1±0.2). Except for the intrinsic z/g 
contamination, which has been deduced from detailed MC simulations, all the background 
components have been measured in different time and energy regimes with the KARMEN 
detector and extrapolated into the evaluation cuts applied for this search. 




RESULTS AND DISCUSSION 

The expected number of z/g induced events from /U+ ^ e+ + z/g + *^z/) is determined by the 
detection efficiencies of the prompt positron and the delayed neutron. The overall detection 
efficiency for positrons is given in Table|l]for a set of different spectral parameters p including 
the contribution from ^^C ( z/g , e^ n ) ^^B, which effectively amounts to less than 5% of the 
p ( z/g , e+ ) n signal in the energy interval of 16-50 MeV. 

Based on the Poisson statistics of the numbers of candidate events and expected back- 



ground, one can extract an energy-independent upper limit for an additional signal [2^ of 
N(z/g)< 7.4 excess events. However, there is additional spectral information, as can be seen 
from Figure Q To use this, we applied a maximum likelihood analysis varying the strength 
of a z^g signal with the energy distribution according to a set of different p parameters. Ta- 
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TABLE I: Flux averaged cross section (a) for p{ue )n and ""^^C n) ^^B, total efficiency 



for e"*" detection, expected {e~^,n) sequences for decaying entirely via /i"* 



e++ + 



experimental results for potential z^g-induced events and deduced upper limits for the branching 
ratio for different spectral parameters p. 
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0.452 
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-0.1 


< 6.2 


< 1.3 • 10-3 


0.5 


86.0 




7.2 


0.456 


5273 ± 489 


-0.4 


< 6.0 


< 1.1 • 10-3 


0.75 


93.5 




8.5 


0.462 


5828 ± 538 


-0.8 


< 5.3 


< 0.9 • 10-3 



blelHshows the signal strength N(z/e)bestfit from the likelihood method. To extract 90%C.L. 
intervals for N(z/e), we performed large samples of MC simulations reproducing experiment- 
like spectra under different signal hypotheses. Our experimental result is consistent with 
no Ue emission from fi'^ decay with upper limits given in Table HI extracted within a unified 
frequentist analysis near the physical boundary N(z/e)=0 following [2^. 

With a potential signal strength of N(z/e)= 5828 ± 538 for a branching ratio BR = 1 for 
LF number violating decays and p(z/e) = 0.75, we set an upper limit of the branching ratio 
of 

> e++ i^e + ^i^^;p{iye 



BR 



0.75) 



e+ + I/e + z/„ 



< 9 ■ 10-^ 



with 90% CL. as well as the upper limits given in TableUfor other values of p(z/e). Figure^ 
shows the visible e^ energies from /i^ — e"*" + z/g + for two different p parameters with 
total strength excluded at 90%C.L.. 

The above limits on the branching ratio BR on /x"*" decays emitting z/g improve by more 
than an order of magnitude the most sensitive limit so far of BR[ii^ — >■ e^ + t'e + < 
0.012 obtained by the E645 experiment at LAMPF 3,Q- 

In models extending the SM, the LF violating muon decay p^ ^ + Ve + '^fi is often 
related to other LF violating processes, e.g. fi —>■ 3e, t pee or muonium-antimuonium 
(MM) conversion. Therefore, limits such as the limit on the probability for spontaneous 
conversion P{M — > M) < 8.2 ■ 10~^^ (90% CL.) j27| set also stringent limits on the coupling 
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constants responsible for the decay /i"*" — ^ e"*" + z/g + 

The most conservative upper hmit of BR < 1.7 ■ 10^^ for any parameter p derived here is 
in direct experimental disagreement with the possibility that the beam excess of Ue seen in 
;he LSND experiment with a branching ratio or probability of P = (2.64 ±0.67 ±0.45) • 10^^ 
lJ | is due to /i^ decays with z/g emission. In the LSND maximum likelihood analysis of the 
data, the best fit to the data is found to be an oscillation contribution with the parameters 
Am^ = 1.2eV^ and sin^(20) = 0.003 In this analysis, oscillation events arise from 

p ( z/g , e"^ ) n with z/g via z/^^ z/g from fi^ DAR as well as from ^^C ( z/g , e~ ) ""^^N with z/g via 
i^fi-^ from vr^ decays in flight (DIF). For large Am^, about 30% of the oscillation signal 
is due to z/^— z/g from DIF, leading to a complex superposition of energy distributions for 
the prompt events distorted by two different LjE combinations, with L being the distance 
source-detector and E the neutrino energy. 

To be able to compare our result quantitatively with the LSND evidence, we follow the 
detailed statistical analysis of the LSND data described in [2^. In this analysis, a special 
cut had been applied to select z/g via z/^— > Ve from yU+ DAR only with almost no oscillation 
events from z/^j— * z/g from vr^ DIF. For large differences of the squared mass eigenvalues 
Am^ = 100 eV^, the extracted interval of the mixing amplitude from the parameter region 
of 90% C.L. was 3.8 ■ lO'^ < sin2(2e) < 8.0 ■ 10-^ with the best fit at sin2(2e) = 5.8 • lO'^ 
corresponding to 65.8(1.3) z/g(z/g) from z/^— > z/g(z/^— > z/g), respectively. For such large values 
of Am^, the energy spectrum of z/g is in good approximation the one given in equ. with 
p = 0.75 2^. Assuming 100% decay probabihty with a spectral parameter p = 0.75, one 
would expect 22692 z/g events seen in the LSND detector. Taking the 65.8 z/g events deduced 
by the maximum likelihood analysis as emitted in the decay /i^ — > e"*" ± z/g ± , we extract 
a branching ratio of BRi^snd{p = 0.75) = 2.9 ■ 10~^ with an approximated 90% C.L. interval 
of 1.9 ■ 10-3 < BRlsnd < 4.0 ■ 10-3. Our derived 90% C.L. limit of BR < 0.9 • 10-^ clearly 
excludes the LSND interval of same confidence and underlines the incompatibility of the 
two experimental results interpretated in terms of LF violating decays. 
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Forschung, the Particle Physics and Astronomy Research Council and the Council for the 
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